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Abstract 

Fragmentation reactions induced on light and medium nuclei by 
protons and light nuclei of energies around 1 GeV/nucleon and be¬ 
low are studied with the Los Alamos transport code MCNP6 and 
with its CEM03.03 and LAQGSM03.03 event generators. GEM and 
LAQGSM assume that intermediate-energy fragmentation reactions 
on light nuclei occur generally in two stages. The first stage is the 
intranuclear cascade (ING), followed by the second, Fermi breakup 
disintegration of light excited residual nuclei produced after the INC. 
GEM and LAQGSM account also for coalescence of light fragments 
(complex particles) up to ^He from energetic nucleons emitted during 
INC. We investigate the validity and performance of MCNP6, GEM, 
and LAQGSM in simulating fragmentation reactions at intermediate 
energies and discuss possible ways of further improving these codes. 


1 Introduction and theoretical background 

Fragmentation reactions induced by protons and light nuclei of energies 
around 1 GeV/nucleon and below on light target nuclei are involved in 
different applications, like cosmic-ray-induced single event upsets (SEU’s), 
radiation protection, and cancer therapy with proton and ion beams, among 
others. It is impossible to measure all nuclear data needed for such applica¬ 
tions; therefore, Monte Carlo transport codes are usually used to simulate 
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impacts associated with fragmentation reactions. It is important that avail¬ 
able transport codes simulate such reactions as well as possible. 

The Los Alamos Monte Carlo transport code MCNP6 [T] uses by default 
the latest version of the cascade-exciton model (CEM) as incorporated in 
its event generator CEM03.03 to simulate fragmentation of light nuclei at 
intermediate energies for reactions induced by nucleons, pions, and photons, 
and the Los Alamos version of the quark-gluon string model (LAQGSM) 
as implemented in the code LAQGSM03.03 (see [2] and references therein) 
to simulate fragmentation reactions induced by nuclei and by particles at 
energies above ~ 3.5 GeV, up to about 1 TeV/nucleon. 

Generally, both CEM and LAQGSM assume that nuclear reactions oc¬ 
cur in three stages. The first stage is the IntraNuclear Cascade (INC), com¬ 
pletely different in CEM and LAQGSM, in which primary particles can be 
re-scattered and produce secondary particles several times prior to absorp¬ 
tion by, or escape from the nucleus. When the cascade stage of a reaction 
is complete, GEM and LAQGSM use the coalescence model to “create” 
high-energy d, t, ^He, and '^He via final-state interactions among emitted 
cascade nucleons, already outside of the target. The subsequent relaxation 
of the nuclear excitation is treated in terms of an improved version of the 
modified exciton model of preequilibrium decay followed by the equilibrium 
evaporation/fission stage of the reaction. But if the residual nuclei after 
the ING have atomic numbers with A < 13, CEM and LAQGSM use the 
Eermi breakup model to calculate their further disintegration instead of us¬ 
ing the preequilibrium and evaporation/fission models. Thus, for targets 
with A < 13, reactions are assumed to occur only in two stages. 

The “standard” version of CEM and LAQGSM account for possible mul¬ 
tiple emission of only n, p, d, t, ^He, and ^He during the preequilibrium stage 
of reactions (see Ref. 12]). Their latest, “E”, version (see Refs. [anE]) 
considers a possibility of preequilibrium emission of light fragments (LE) 
heavier than ^He, up to ^®Mg. It also simulates coalescence of LE heavier 
than ^He, up to A = 7, in CEM03.03E (see [3l[ll[5]), and up to A = 12, in 
LAQGSM03.03F (see [6|). 

In recent years, MCNP6, with its CEM and LAQGSM event generators, 
has been extensively validated and verified (V&V) against a large variety of 
nuclear reactions on both thin and thick targets (see, e.g. Refs. lal-lEl and 
references therein). In Ref. [3|, it was tested specifically on fragmentation 
of light nuclei at intermediate energies. Here, we present a few results from 
our recent work [3| and investigate further the performance of MCNP6, 
CEM, and LAQGSM in simulating fragmentation reactions at intermediate 
energies and discuss possible ways of further improving these codes. 


2 Results and conclusion 


Figs. 1, 2, 3, and 5 show examples of fragmentation reactions on light nuclei 
simulated by our codes. Figs. 6 and 4 shows examples for medium targets, 
^®Ca and ”“*Ag. Many more similar results, their discussion, and useful 
details can be found in Refs. m-M and references therein. 

Figure 1: Excitation function for the 
production of from p + calcu¬ 
lated with CEM03.03 using the “stan¬ 
dard” version of the Fermi breakup 
model {Apermi = 12) and with cut¬ 
off values for Apermi of 16 and 14, as 
well as with MCNP6 using CEM03.03 
{Apermi = 12) Compared with experi¬ 
mental data, as indicated. Experimen¬ 
tal data are from the T16 Lib compila¬ 
tion [To] (see details in B)- 


Figure 2: Comparison of experimental 
®Li spectra at 20, 45, 60, 90, and 110 
degrees by Machner et al. [TT] (sym¬ 
bols) with calculations by the unmod¬ 
ified CEM03.03 (dashed histograms) 
and results by CEM03.03.F (solid his¬ 
tograms) , as indicated (see more details 
in [3]). 

Our results indicate that 
MCNP6 using CEM03.03 and 
LAQGSM03.03 simulates fragmen¬ 
tation reactions on light and 
medium-light nuclei at intermediate energies well, in a satisfactory agree¬ 
ment with experimental data. The recent “F” version of codes (see Refs. 
[3] - [6]) is even better, as it allows us to describe emission of energetic LF 
from practically arbitrary reactions. 

However, MCNP6 is not yet ready to predict well heavy fragments from 
reactions with heavier nuclei, with mass numbers A ~ 100. Such nuclear tar¬ 
gets are considered too light to fission in conventional codes. Similarly, the 
fragments are too light to be produced as spallation residues and too heavy 
to be produced via standard evaporation and/or preequilibrium models, or 
via coalescence. 
























Figure 3: Comparison of measured 

[H] t spectra at 45, 60, 90, and 130 
degrees from 800 MeV/nucleon ^°Ne + 
^°Ne with calculations by LAQGS03.03 
using its “standard” version of the co¬ 
alescence model {po = 0.108 GeV/c 
for t and ^He; dotted lines) and with 
modified values of po labeled in legend 
as “coall” (dashed lines) and “coal2” 
(solid lines), as indicated in legend and 
discussed in detail in Ref. [3]- 


Figure 4: Comparison of experimen¬ 
tal data by Green et al. „II3] ( cir¬ 
cles) for the production of ®Li at an 
angle of 60° from the reaction 480 
MeV p -I- "-“‘Ag, with results by stan¬ 
dard CEM03.03 (dot-dashed line) from 
CEM03.03F without coalescence ex¬ 
pansion (solid line) and CEM03.03F 
with coalescence expansion (dashed 
line) (see details in [5lfT^l. 

One way to approach this problem 
would be to employ after the INC 
stage of reactions a fission-like 
sequential-binary-decay model, like the code GEMINI by Charity et al. m 
to describe the compound nuclear decay. In our case, this means sep¬ 
arately merging GEM and LAQGSM with GEMINI. Actually, we already 
have done so more than a decade ago, producing the “G” versions of GEM 
and LAQGSM we had at that time (see, e.g., Ref. m and references 
therein). 

Another way to address this problem is to implement in GEM and 
LAQGSM the Statistical Multifragmentation Model (SMM) by Botvina et 
al. m- Thus, we would consider multifragmentation as a mode competitive 
to evaporation of particles and light fragments, when the excitation energy 
E* of a compound nucleus produced after the preequilibrium stage of a re¬ 
action is above a certain value, e.g., = 2 x A MeV, as we did in 

the “S” versions of CEM03.01 and LAQGSM03.01 (see, e.g.. Ref. [E] and 
references therein). 

As of today, neither the “S” nor the “G” versions of CEM and LAQGSM 
have been implemented in MCNP6. We plan to incorporate them in our 























event generators used by MCNP6 after we tune several parameters in SMM 

Figure 5: '*He spectra (at 35°) from 
1.2/1.9/2.5 GeV p + measured 



to multifragmentation simulated with 


by M. Fidelus of the PISA collabora¬ 
tion [13] (symbols) with calculations by 
MCNP6 using CEM03.03 (see details in 

0 ). 

and GEMINI, that are essential 
in chosing the excitation energy 
(or temperature) of nuclei when 
reaction mechanisms change from 
“usual evaporation”, to binary de¬ 
cays described by GEMINI, and/or 
SMM (see details in m)- 


Figure 6: Measured cross sec¬ 
tions for ^®Ca fragmentation on 
®Be at 140 MeV/nucleon [T5] 
compared with LAQGSM03.03 
predictions (see details in [9]). 
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